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Recent experiments showed that the surface of a three dimensional topological insulator develops 
gaps in the Floquet-Bloch band spectrum when illuminated with a circularly polarized laser. These 
Floquet-Bloch bands are characterized by non-trivial Chern numbers which only depend on the 
helicity of the polarization of the radiation field. Here we propose a setup consisting of a pair of 
counter-rotating lasers, and show that one-dimensional chiral states emerge at the interface between 
the two lasers. These interface states turn out to be spin-polarized and may trigger interesting 
applications in the field of optoelectronics and spintronics. 


PACS numbers: 73.20.At; 78.67.-n; 73.43.-f; 72.25.-b 

Introduction. - Amid the thrill sparked by graphene [TJ 
[ 2 ] and its record properties [3], the discovery of topo¬ 
logical insulators (TIs) [H [5] developed with surprising 
speed. Indeed, TIs were predicted two years earlier in 
graphene [6], but the necessary spin-orbit interactions 
were too weak for this to be observed and a different 
playground was needed to realize them [ 7 ] . Most TIs are 
three-dimensional materials like usual solids, but with a 
special property: they have a bulk band gap while keep¬ 
ing states that propagate with unprecedented robustness 
at the periphery of the sample 0 0. These peculiar 
states hold great promise for quantum computation m 
but at the same time open up a major challenge: con¬ 
trolling them is particularly demanding for 3D TIs. 

Encompassing the rapid progress in graphene photon¬ 
ics m and optoelectronics [12] [13], theoretical studies 
predicted the formation of laser-induced band gaps m 
in graphene when properly tuning the laser polariza¬ 
tion, frequency and intensity [TbHTS] . More recently, 
these gaps were unveiled at the surface of a TI through 
ARPES [19]. This triggered great expectations for 
achieving laser-assisted control not only in the form of 
an on-off switch for the available states but also because 
theoretically non-trivial topological states [Huniiii] can 
be induced on a diversity of materials [22H26] . and also 
in cold matter physics [27] [28] . Exciting questions arise 
about the nature of these novel states [2911^ , the possi¬ 
bilities for manipulating them m, the associated topo¬ 
logical invariants [32H36] , their statistical properties [S3- 
SO] and their two-terminal gU HU and multi-terminal 
(Hall) response [43l [44]. Still, an experimental realiza¬ 
tion of the Floquet chiral edge states is missing. Most 
studies considered two-dimensional systems, except for 
Refs. [45] and [46] where the target was a 3D semicon¬ 
ductor. 

Here we show that besides opening a band gap as in 
Ref. [19] , illuminating a 3D TI with a suitable set of lasers 
can confine the surface states into one-dimensional states 
which also bear a topological origin. The proposed setup 



FIG. 1. (color online). Scheme of the proposed setup where 
the surface of a 3D TI is illuminated by two circularly polar¬ 
ized lasers with opposite directions of rotation. The y-axis de¬ 
fines the interface region between the two lasers. The dashed 
line in the border indicates that the system extends indefi¬ 
nitely in the x-y plane. Inset: Laser intensities as function 
of the x-coordinate. At the interface, the two lasers interfere 
and the total field becomes linearly polarized. 

is represented in Fig. Two lasers with opposite circu¬ 
lar polarization incident perpendicularly to a face of a 
3D TI. This can be devised through, e.g., a single laser 
beam splitted into two with opposite helicity. The in¬ 
terface between the lasers is assumed to be shorter than 
the thermalization length so that the occupations are de¬ 
termined by the (larger) regions without radiation. As 
we will see below, this modification of the experimental 
setup in Ref. [19] introduces Floquet states propagat¬ 
ing along the boundary where the polarization changes. 
Our results follow from simulations of the Floquet spectra 
based on low-energy models, which are further supported 
by: (z) a calculation of the topological invariants and {ii) 
explicit calculations for a driven 3D lattice model. Inter¬ 
estingly, we show that the resulting Floquet boundary 
states, which arise from a topological transition between 
the illuminated regions, carry spin-polarized currents. 

Illuminated TIs and Floquet theory.- We consider a 
low-energy Hamiltonian describing the surface of a TI. 
Assuming the (001) direction [47] and linear order in k, 
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the effective surface Hamiltonian reads Hq = hv{ky(7x — 
where cTx and cFy are Pauli matrices describing 
the spin degree of freedom. The time-dependent field 
is included through the Peierls substitution k ^ k ^ 
eA{t)/hc^ with A{t) the laser’s vector potential. In the 
regions dominated by one of the two lasers, i.e. \x\ ^ xq, 
with xq the characteristic length of the lasers’ interface 
(see Fig.[^, we choose a circularly polarized field Ar{t) = 
Aq [cos(rOt ex-\- sin(rOt + cp) By] , where r = ±1 sets 
the direction of rotation and cp determines its orientation 
(measured from the x-axis) at t = 0 and, as shown later 
on, it becomes relevant at the interface’s region x ^ 0 . 
The time-dependent Hamiltonian thus reads 

Hr{t) = Hq + 7 sin(r(2t + (/:?) ax — j cos{Tflt ^ (p) ay, (1) 

where 7 = vcAq/c characterizes the strength of the per¬ 
turbation. A suitable description of the dc-spectrum and 
the topological properties of the system can be achieved 
through the Floquet theory. By using Floquet’s theo¬ 
rem, we obtain a time-independent Hamiltonian in Flo¬ 
quet space, defined as the direct product IZ^T between 
the usual Hilbert space IZ and the space of time-periodic 
functions T. This space is spanned by the states |Tcr, m), 
where a = {t, 1} accounts for spin and m is the Fourier 
index. The Floquet Hamiltonian writes 

n],ik) = Ho0l+imn+nT ( 2 ) 

/3=± 

where we use a± = {ax ± iay)j2. Such Hamiltonian 
can be imagined as a series of replicas (Floquet chan¬ 
nels) of i^o, each one defined in a Fourier component 
of the driving. The static i^o enters in the diagonal 
part, together with the contribution [A^o]n,m = mhClSn,m 
from the driving field, and the vector potential couples, 
through [A/ 3 ]n,m = thosc channels differing in 

their Fourier indices by Am = ± 1 . 

For the calculation of the laser-induced band gaps and 
the associated Chern numbers, it is enough to consider an 
homogeneous system defined at the TI’s surface through 
Eq. Q. The underlying assumption is that HQ is smaller 
than the bulk gap, such that the states associated to the 
bulk do not participate in the gap openings. As discussed 
in Refs. [22l|35], these laser-induced gaps are indeed de¬ 
pletions of the time-averaged density of states which re¬ 
sults from weighting the Floquet spectrum on the m = 0 
channel. By assuming low intensities {j/hfl 1) we re¬ 
strict ourselves to the main contributions to the band gap 
openings around 5 = /iO /2 and at the Dirac point 5 = 0, 
henceforth called the zone boundary (ZB) and the zone 
center (ZC) gaps, respectively. These two gaps were de¬ 
scribed in Ref. [48], obtaining Ai ^ j and Aq ~ jHQ 
for the ZB and ZC gaps, respectively. Notice that a 7 r/ 2 - 
rotation along the 2 ;-direction of the spin coordinate sys¬ 
tem maps Hq to the low-energy Hamiltonian describing 
a single valley in graphene {Hq hvk • cr). Therefore, 


apart from a change in the Fermi velocity, the laser- 
induced gaps show the same dependencies in both sys¬ 
tems [11 [13 Eg. 

The equivalence between Eq. (§ and the low-energy 
description for illuminated graphene can be exploited 
even further: In graphene, the laser-induced gaps are 
characterized by non-trivial Chern numbers, and the 
bulk-boundary correspondence leads to Floquet chiral 
states at the sample edges [H ESI ED- Can similar states 
appear here? A first problem is simply that the surface 
of a 3D solid cannot have a boundary. This motivates 
our proposal of changing the light polarization as in Fig. 

thereby introducing an effective boundary (much like 
a domain wall, as discussed in other examples [5Qll52] i 
where Floquet chiral states develop—by chiral we mean 
that the direction of motion is fixed by the helicity of the 
two lasers. 

Starting from Eq. Q, we proceed as in Refs. [351153]: 
Eirst we isolate each crossing where a band gap opens, 
and then we compute a 2 x 2 effective Hamiltonian of 
the form = h cr. The contribution to the Chern 
number of the lower band can be calculated through the 
expression [9]: 

( 3 ) 

In the ZB gap region, the band gap opening comes from 
the crossing between the states |T+, 0) and |T_, 1). Here 
|T±) refer to the conduction and valence band solutions 
of Hq, respectively, and the second index (0 or 1) in¬ 
dicates the Eloquet channel. By reducing the Eloquet 
Hamiltonian to these states, we obtain that the contri¬ 
bution to the Chern number is Ci = r. In the ZC re¬ 
gion there are two processes related to the gap opening. 
These consist of (z) the renormalization of |T±,0) due 
to the coupling to the m = ±1 states, and {ii) the level 
crossing between |T_,1) and |T+,—1), bridged by the 
m = 0 states. A straightforward calculation of these two 
contributions yields Cq = — t/2 + 2t = 3r/2. While in 
graphene this half-integer Chern number is compensated 
by spin and valley degeneracies, in strong TIs, where the 
surface encloses an odd number of Kramers degenerate 
Dirac points, a half-integer Chern number results for ex¬ 
ample when exposing the material to a static magnetic 
field [illlllllHk]. 

Interface states in 3D TIs.- A natural question relies 
on the bulk-boundary correspondence in illuminated TIs 
associated to the non-zero Chern numbers of the Elo¬ 
quet bands. In the present case, since inverting the he¬ 
licity of the circularly polarized laser changes the sign of 
the Chern numbers, one expects the generation of chiral 
states at the boundary between the two regions. 

To elucidate this question we proceed by solving the 
proposed model at the laser’s interface. Eor simplic¬ 
ity in the calculation we assume a sudden change of 
the laser’s direction of rotation by assigning a different 
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FIG. 2. (color online). Laser generation of interface states 
crossing the ZB (a) and ZC (b) gaps. Gray regions indicate 
extended states zones. The color scale shows the time-average 
spin texture (((7a;)). Inset: Time-averaged probability density 
P for the state at £ = HQ/2 in solid line (black). Dotted 
line (green) shows the m = 0 component while dashed line 
(purple) is for m = 1. Here we use — 0.2 and if — 7r/2. 

r to each portion of the system [according to Fig. 
t{x) = —sgn(x)]. The resulting differential equation 
therefore reads 

d^^{r) = Mr(x)^{r), (4) 

where Air = ^y) ~ and ^(r) = 

contains the wave-function 
coefficients for the involved channels [58]. The 

solutions to the above differential equation follow a stan¬ 
dard diagonalization of Afr in the two regions with the 
appropriate boundary condition m and is discussed in 
detail in the Supplemental Material [59] . 

In Fig. we show the resulting localized states in both 
the ZB and ZC gaps in a configuration of f = 7r/2, 
such that at t = 0 the vector potentials point parallel 
to the interface’s direction. Our calculations allow to 
verify the bulk-boundary correspondence. Indeed, the 
difference between forward and backward propagating 
states (relative to the ^-axis) is always fixed to the cal¬ 
culated difference between the Chern invariants at each 
side of the interface, yielding ACi = —2 for the ZB gap 
and ACq = —3 for the ZC gap. Other choices in the 
orientation cp lead to changes in the dispersions of the 
chiral states but keeping these numbers the same m- 
Although never crossing the gaps, additional states lo¬ 
calized at the interface are also present—since they are 
not chiral, one could expect these states to backscatter 
when encountering impurities in the sample. In Fig. ia) 
(inset) we show the time-averaged probability density. 


P{x) = a associatcd to one of the two lo¬ 

calized states crossing the gap at e = hft/2 (negative ky). 
Here, the decay length depends inversely on the size of 
the gap and, similarly to illuminated graphene [221135], 
it is independent of the microscopic details of the sam¬ 
ple. In a more realistic situation where the inversion of 
the laser helicity is taken gradually over a finite length 
xq (see below), the width of these states shows to de¬ 
pend also on the latter parameter xq. It can be seen 
also that there is a pronounced asymmetry in the con¬ 
tributions from the m = 0 and m = 1 channels to the 
overall probability density, which is particular to the rel¬ 
ative angle between p and the direction of the interface. 
The other state developed at positive ky shows to have 
this asymmetry inverted. 

Using the solutions of Eq. Q we calculate the time- 
averaged spin texture [35] associated to the Floquet 
boundary states. Thanks to the spin-momentum lock¬ 
ing present in the TI without radiation, there is a non¬ 
vanishing spin component in the Floquet states, {{<Jx)) = 
2 / Re['0|^^(x)'0^,^(x)]dx, i.e. the in-plane compo¬ 

nent perpendicular to the interface’s direction. In all 
cases, {{(Jx)) is proportional to the group velocity, as can 
be seen in Fig. Since the direction of propagation can 
be tuned by the lasers, these robust states could be also 
interesting from the point of view of spin-polarized trans¬ 
port at a desired region of the TI’s surface. 

Three-dimensional lattiee model and LDOS.- Up to 
now our analysis is based on an effective 2D-model for the 
surface states of a TI. This poses the question on whether 
these properties can be reproduced in a model accounting 
for the insulating bulk bands of a 3D TI. We therefore 
consider a lattice Hamiltonian which satisfies the four 
symmetries present in a strong TI [47]. By taking a cu¬ 
bic lattice with parameter a, we obtain a tight-binding 
description for an isotropic TI m- The vector potential 
A{r^t) = Ar{z,t)f{rx)^ enters through Peierls’ 

substitution as a time-dependent modulation of the hop¬ 
ping matrices coupling nearest neighbor sites. This ac¬ 
counts for (i) a gradual change f{x) = [1 +exp(x/xo)]“^ 
of the laser helicity, which produces a (^-oriented, lin¬ 
early polarized field at the interface region, and {ii) a 
photon absorption process across the layers of the TI 
which manifests through a decay in the laser intensity, 
Ao{z) = In onr simulations, ^ and zq are ad¬ 

justed in such a way that the laser becomes negligible 
at the bottom face of the irradiated sample. The result¬ 
ing lattice Hamiltonian is derived in the Supplemental 
Material [59] . 

In Fig. we show the time-averaged /c^-resolved LDOS 
(pky) [mlS] in a geometry [see panel (a)] where the solid 
is infinite along the x and y directions, while in the 2 :- 
direction it has Nz = 9 layers. A quantitative description 
is possible by adjusting the model parameters to, e.g., 
those estimated in Ref. [60] , yielding larger values of Nz . 
In panel (b) we calculate pky for the non-illuminated ma- 
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FIG. 3. (color online). Normalized pky in an isotropic model for a 3D TI with Lz = 9a: (a) Schematics of the considered 
setup indicating the different points at which pky is evaluated: (b) At the surface of a static TI. (c) Irradiated sample at the 
insulator’s bulk {z = 5a). (d) At the surface, away from the interface’s region, here only one of the two lasers dominates, (e) 
At the surface, across the interface’s region, (f) Used color scale in all plots and close-up of panel (e) at the ZB (top panel) 
and the ZC (bottom panels) regions. In all cases, the used parameters for the static Hamiltonian are mo — 0.4, mi = 0.25, 
and m 2 = 0.5 [51], and the laser parameters are hQ = O. 270 , po = 27raoa/4>o = 0.1 [59], p = 7r/2, xo = 3a, ^ = 0.8, zo = a and 
includes the Floquet channels ?n = { —1,0,1}. 


terial, where the gapless surface state crossing the bulk 
gap can be appreciated. Turning on the lasers, we evalu¬ 
ate pky at different points of the sample. Panel (c) shows 
the sample’s bulk region, where there is a clear absence 
of states within the insulating bulk gap. When mov¬ 
ing to the top surface, to a region dominated by only 
one of the two lasers, panel (d) reveals the ZB and ZC 
gaps similar to those observed in Ref. m- Finally, once 
we arrive to the interface region, panel (e) shows the 
emergence of chiral states similar to those of Fig. In 
the ZC region [bottom panels in Fig. [^f)], we can ob¬ 
serve that due to the small width of the gap the central 
state (forward mover) apparently crosses it, reflecting the 
—t/2 contribution to Cq from the renormalization of the 
m = 0 states. Similar to Fig. [^b) there is, however, 
a small admixture which hybridizes this forward mover 
state with the degenerate states around ky ^ 0 (back¬ 
ward movers) and the final state crossing the gap shows 
a negative slope, as required by ACq. In this sense, the 
difference in the number of states with opposite direction 
of motion is again bounded to the calculated topological 
invariants Co and Ci on each side of the interface and do 
not depend on its specific shape. The details of the wave 


functions and of the quasi-energy dispersion, however, do 
depend on the angle between the interface and the orien¬ 
tation p of the linearly polarized vector potential formed 
at that point m- 

Final remarks. - In summary, we found that illuminat¬ 
ing the surface of a 3D TI with a pair of counter-rotating 
lasers generate chiral, one-dimensional states confined 
at the interface region between the lasers. These states 
locate within the recently measured laser-induced gaps 
in ARPES [19], for which we believe a small modi¬ 
fication of the experimental setup would be enough 
for its observation. Additionally, these states have 
a finite time-averaged spin texture subjected to the 
spin-momentum locking effect of the bare material, 
making them interesting from the point of view of spin 
polarized transport. Our calculations in the low-energy 
regime are supported by simulations in a 3D lattice 
model, which accounts for the interface zone. Given 
the topological character of the Floquet bands, the 
qualitative properties of these interface states (chirality 
and spin-momentum locking) remain unaffected by the 
experimental details of the laser configuration {e.g. 
fluctuations in their relative phase). Importantly, the 
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existence of the topological states is not bounded to 
the local recovery of the time-reversal symmetry at the 
interface m- Other choices in the setup including, e.g.^ 
the simultaneous irradiation of different faces of the TI, 
are of great interest and deserve further exploration 
since one could exploit different spin-textures and band 
curvatures ga EH to achieve control over the chiral 
states. 
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SUPPLEMENTAL MATERIAL 

In this supplemental material we provide additional 
details concerning: (a) The solutions of the differential 
equation (continuum model) presented in the main 
article, (b) the explicit form of the Floquet Hamiltonian 
for the lattice model of an isotropic 3D TI, (c) the role 
of the linear polarization angle ip in the dispersion of the 
laser-induced interface states, and (d) the case of two 
lasers with different frequencies. 

Solutions of the continuum model.— The discussed 
differential equation in the main article for an interface 
region along the ^-coordinate is of the form 

a,^(r) = 7W,(,)^(r), (5) 

where in general ^(r) = (...,with 

'04,,m)^: accounts for spin-up and spin-down 
states for each one of the m-Floquet channels and the co¬ 
efficient matrix Afr = ~ with de¬ 
fined in Eq. (§, comes from the replacement —idx 

due to the broken traslational invariance. The solutions 
of Eq. <§ can be obtained by diagonalizing m. the 
two regions {x ^ 0) of the sample separately. These are 
of the form 


Ux) = (6) 

3 

where i labels the states of the truncated basis {a, m}, 
Ur is the transformation matrix that diagonalizes Mr 
and t{x) = —sgn(x) determines the direction of rotation 
of the vector potentials at each side of the interface. For 
each one of the considered gapped regions, namely the 
zone center (ZC) and zone boundary (ZB), we work in a 
different truncation basis of Floquet channels to ensure 
a symmetric eigenvalue spectrum of Mr around zero. 


Specifically, we work in the Floquet space defined by the 
m G {0,1} for the ZB gap and m G {—1, 0,1} for the ZG 
gap. This last guarantees that for an eigenvalue \j of Mr 
there is always another A/, such that = —Xj and allows 
us to order them in the form Re(Ai )<...< Re(A7v/2) < 
0 < Re(A7v/2+i) < ••• < Re(AAr), with N the dimension 
of the truncated space. Due to the specific form of the 
coefficient matrices Mr, the eigenvalues Xj are indepen¬ 
dent of the r-index, yielding the same A-spectrum in the 
two regions [62]. To ensure the convergence of for 
X Too in the considered gaps, we set to zero those 
coefficients which are associated to Re(Aj) ^ 0. Ac¬ 
cording to the above ordering, this implies 

C+ = 0, j = l,...,7V/2, (7) 

C-=0, j = 7V/2 + l,...,7V. (8) 

The remaining coefficients are found by imposing a topo¬ 
logical boundary condition m across the interface be¬ 
tween the two portions of the system. In the considered 
setup, this implies the continuity of ^(r) along the x- 
direction, where the sign in the Chern number is inverted, 
and it reads 


N 

^ [ur^ 


i=JV/2+l 


+ \^3 


N/2 

c; = 

i=i 


ij^j 


(9) 


Since in the above equation there is only one coeffi¬ 
cient for each particular state j, we can define C = 
{C^ ,...,CM/ 2 XN/ 2 +v-XtiP, such that QC = 0 , with 


roi ^ r -[wrE. i = 

1 +[Up]ir j = 7V/2 + l,...,7V. 


( 10 ) 


Through the condition det Q = 0 we thus determine 
numerically the energy 5 and momentum ky of the 
interface states within the ZC and ZB gaps where all the 
A-eigenvalues have a non-vanishing real component (cf. 
Figs. and 1^. 


Spin texture.- Based on the above solutions for the 
interface state wave-functions we can now calculate the 
expectation values of the spin operator a. Let’s assume 
we find a solution ^(r) such that it satisfies the charac¬ 
teristic equation det Q = 0. Written in the Hilbert’s real 
space (7^), the wave-function reads [35] 

\^{r, t)) = E , (11) 

a,m 


where {|Ter)} is a complete basis for the spin states in IZ 
and the expectation value of the spin operator therefore 
reads 

(o-) = E E 

cr,cr' m,m' 


For time-averaged quantities over a period T = 27r/(7 of 
the driving we observe that only the direct terms with 
m' = m survive and hence 

((-)) = EE/ 

a,a' m 

This last was calculated in Figs. and [^for the interface 
oriented along the ^-direction and yields the same 
spin-momentum locking effect observed in static 3D TIs, 
i.e., {{(Ty)) = ((cr^)) = 0 and finite {{<Tx)). 


Lattice Floquet Hamiltonian.— Here we derive the 
Floquet Hamiltonian for the lattice model introduced in 
the main article. We start from the static Hamiltonian 
for the isotropic TI of Ref. m for a cubic geometry with 
lattice parameter a: 

^ = E clMoCr + E 4r„Cr+oe„ + H.C., (13) 

r r,a 


where the sum runs over the lattice’s sites r and a = 
{x, z}. The 4x4 on-site Mq and hopping matrices 

Mo = (mo - 6mi)T^, = (miT;, - , (14) 

account for both the orbital and spin degrees of free¬ 
dom through the Pauli’s matrices and (7^, respectively. 
Here mo, mi and m2 are standard parameters of the 
model [51] whose scale is set by the hopping term 70. 
The driving field is included by Peierls’ substitution as 
a phase-modulation of the hopping amplitudes coupling 
nearest-neigbour sites. This incorporates space and time 
dependencies in through the vector potential, defined 
as A(r,t) = — /(x)), where 

f{x) = 1/[1 + exp(x/xo)] and 


At- = Ao{z) [cos(rr^t -j- cp) Cx sin(rOt + cp) Cy] , (15) 

with Ao{z) = the magnitude of the vector po¬ 

tential, assumed to be attenuated along z due to absorp¬ 
tion within the solid. At the interface region {x = 0) 
the two lasers add-up, yielding a linearly polarized field 
A = Ao( 2)) cos(flt)e^ whose direction forms an angle ip 
with respect to the x-axis. The phase-modulation thus 
enters through the line integrals 


Ta{r,t) -)■ T^exp 


2 ^ rr^aeoc 

i— / A{r,t)-d£o 

Jr 


, (16) 


with 4>o the magnetic flux quantum. Now that we have 
the explicit time-dependence of the hopping matrices, the 
Floquet Hamiltonian can be easily obtained by taking the 
Fourier components of the above phase-modulation, i.e., 
[^F]n,m = ^ Jq such that 

'Hf = y]4A4oCr + y^^clTa{r)cr+aea, + H.c. (17) 
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The above static matrices Mq and generalize in Flo¬ 
quet space to 

Mo = MqiSiI+ I®Nq, Ta{r) = Ta®'ya{r), (18) 

where [A^n]n,m = nMl5n,m- The tensor product fol¬ 
lows from the TZ0T structure of the Floquet space and 
[la{r)]n,m = 72 “’"(^)> with 

7^ = E J„(^5) J„_„(r,(l - g)), (19) 

m 

7 = E - 4), (20) 

m 

and 7^ = Here, is the n-th order Bessel 

function of the first kind, g{x) = ^ f{x')dx' and 
r]{z) = 27raAo(^)/4>o, such that 70 = 27raoa/4>o sets the 
strength of the couplings between the Floquet channels. 

Polarization angle.— In this section we describe in 
detail the role of the angle ip of the linearly polarized 
laser formed in the vicinity of the interface region [see 
Figs. and l^a)]. 

We start our description based on the solutions dis¬ 
cussed above for the low-energy, continuum model. In 
this case, we assumed a sudden change in the helicity 
of the laser through r{x) = —sgn(x) [see Eq. (|^], and 
the relevant angular dependence comes from the sum 
between the two contributions, i.e., ^+(^) + ^-(^) = 
Aq cos{Qt)ei^, which yields a linearly polarized field 
whose oscillation direction forms an angle ip with respect 
to the x-axis. Although this linearly polarized field is 
not explicitly present in the continuum model, the dis¬ 
persions of the interface states strongly depend on the 
relative orientation of the two fields which is kept fixed 
by ip at any value of time. 

In Fig. we show the solutions to det Q = 0 discussed 
above in the vicinity of the ZC and ZB gaps. In panel 
(b) we fix the energy at the center of each one of the 
gaps and evaluate the /c^-position of the interface states 
for arbitrary orientation ip. It can be seen that for the 
ZB region (top panel) these states are 7r-periodic while 
in the ZC region (bottom panel) they turn to be 7r/2- 
periodic. As we will discuss below, the difference in the 
(^-periodicity of the two regions bears a direct resem¬ 
blance with the Floquet spectrum in a sample irradiated 
with a linearly polarized field. In panels (c) and (d) we 
explore the dispersion relations of the interface states for 
intermediate angles, where the dependence with ip be¬ 
comes evident. For the ZB gap, cf. Fig. Qc), the two 
interface states move backwards with respect to the in¬ 
terface’s direction y. These start at maximum separation 
for ip = 7r/2 (left panel) and then they shift to the cen¬ 
ter ky ^ d, becoming almost degenerate ioi ip = d (right 
panel). In this (/:?-evolution we observe, in addition, a 
change in the group velocities, and together with it a 
decreasing spin polarization. For the ZC gap [see panel 
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FIG. 4. (color online). Interface states in the ZB and ZC gaps: (a) Schematics of the laser’s rotation along the a:-direction. 
In the vicinity of the interface region (white) the resulting linearly polarized field points in a direction ip with respect to the 
x-axis. (b) Angular dependence of the interface states evaluated at the center of the ZB gap {e = hQ/2, top panel) and ZC 
gap (e = 0, bottom panel), (c) ZB interface states evolution for (p — 2 (left), p — ti/A (center), and p — ^ (right), (d) 

ZC interface states evolution for p = 7r/2 (left), p = Stt/S (center), and (/? = 7r/4 (right). In (c) and (d) the insets (arrows 
and circles) show the angle p and the color scales indicate the time-average x-component of spin {{(Jx)) [cf. Eq. ( |12[ )]. As in 
Fig. 1^ of the main text, we use ^/hTL — 0.2. 


(d)], we notice a central state moving forwards and dom¬ 
inated by a negative spin component. This state remains 
almost in the same /c^-position regardless the value of p. 
The other states (backward movers) shift in such a way 
that those in the center are degenerate for (^ = 7r/2 (left 
panel) and then they merge with those initially placed 
in the extremes for p = 7r/4 (right panel). This marked 
difference between forward and backward movers has in¬ 
deed its origin in the different contributions to the Chern 
number Cq = — t/2 + 2r discussed in the main article. 
These contributions are not, however, completely ‘sepa¬ 
rated’ when including the interface boundary, manifest¬ 
ing itself as an avoided crossing between forward and 
backward movers. In any case, this admixture between 
states with different directions of propagation is in agree¬ 
ment with the bulk-boundary correspondence, since this 
last dictates not the total number of interface states but 
the difference between forward and backward movers. 

To elucidate the /c^-position of the interface states 
in the above discussed dispersion relations, we show in 
Fig. i a) the calculated Floquet spectrum of a sample 
which is being irradiated by a single laser with linear 
polarization [48]. In this setup both the traslational in¬ 


variance along X and the time-reversal symmetry are re¬ 
covered. The corresponding Floquet Hamiltonian in this 
case writes: 

Hf = i7oC)/ + /(8)A^Q + ^(sin(/:?cra,-cos(^cr^)(8)A, (21) 

where Hq = hv{ky(Jx — kx<7y)^ 7 = veAojc and = 

From Fig. j^a) it can be seen that the 
aforementioned laser-induced gaps now close at certain 
points as a consequence of the reinstated time-reversal 
symmetry. This ‘closing’ of the gaps occurs at two points 
in the ZB region (top) while in the ZC region (bottom) 
we have a central Dirac cone (gapless) surrounded by four 
closing points in a different shell of the spectrum. From 
Eq. (21) we notice that the role of the polarization angle 
p is to rotate the whole spectrum and with it the points 
where the gaps close. Returning to the initial setup with 
the two counter-rotating, circularly polarized lasers, we 
can now associate the /c^-position of the interface states 
bridging the gaps to these points in which the gaps close 
when accounting for the linear polarization scenario. In 
particular, in the ZC region the difference between for¬ 
ward and backward movers can be attributed to the shells 


































ky {n/a) 
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FIG. 5. (color online), (a) Typical dispersion relations for the irradiated system with linearly polarized held around the ZB 
(top) and ZC (bottom) regions, (b)-(c) Normalized pky for the isotropic lattice model in Eq. (|17|): (b) In the vicinity of the 
ZB gap, for p — 7r/2 (top), (/? = 7r/4 (center), and p — ^ (bottom), (b) pky around the ZC gapfor p — 'nj 2 (top), p — Z'k 
(center), and p — 'kIA (bottom). The remaining parameters coincide with those used in Fig. of the main article: Nz — 9, 
rriQ = 0.4, mi = 0.25, m2 = 0.5, HQ = O.270, ?7o = 27raoa/4>o = 0.1, xq = 3a, ^ = 0.8 and zo = a. 


of the spectrum to which they belong: The forward mover 
is tied-up to the central Dirac cone, related to the m = 0 
channel, while the backward movers are placed in the 
outer closing points, associated to the crossing between 
m = — 1 and m = 1 cones. Within this correspondence 
between interface states and the points where the gaps 
close, it is possible to understand the above periodici¬ 
ties in both the ZB and ZC gaps [see Fig. [^b)] as a 
consequence of the rigid rotation of the spectrum when 
sweeping p. 

To further support the above discussed (/:?-dependence 
in the dispersion relations of the interface states, we also 
present calculations of the time-averaged LDOS (pky) 
evaluated along the interface region for the lattice model 
described in the previous section. In Figs. ib) and 
(c) we show the normalized pky in the ZB and ZC gap 
regions for the same angles p as in Figs, [^c) and (d), 
respectively. All the discussed effects are reproduced 
qualitatively. In this case where the intensity of the laser 
is smaller than the one used in the continuum model, the 
mixing between forward and backward movers becomes 
almost unnoticeable. 

Lasers with different frequencies.— As a concrete 
test example to inspect the robustness of the laser- 
induced interface states against possible breaking of the 
local time-reversal symmetry (TRS) we consider the case 
of two laser beams defined with two different frequencies 
and opposite helicities. To ensure the global periodicity 
of the radiation field, we take commensurable frequencies. 


Indeed, the TRS is broken everywhere in this system. 

In Fig. I^a) we calculate the determinant of the Q- 
matrix as defined in Eq. ( [1Q| ) around the Floquet zone 
center (ZC) region. The solutions to Eq. (§ correspond 
to those points (black traces) where the map goes to zero. 
Here, we can appreciate the presence and chirality of the 
laser-induced interface states for the case of two laser 
beams with frequencies Qi = Q and O2 = 20, respec¬ 
tively, and opposite helicities. We have tested other cases 
with smaller frequency ratios, and they all show the same 
features but become technically demanding [63]. The 
spectrum results similar to that of the single frequency 
case (cf. Fig. Hb) in the main manuscript). At very 
low energies, differences are expected from this situation 
regarding the effect of additional Fourier replicas that 
need to be included in our calculation. These, however, 
are well-understood and follow the hierarchy recently in¬ 
troduced in Ref. [53] . 

We emphasize that the formation of chiral interface 
states even in this system with fully broken TRS can 
be explained by the closing of the laser-induced ZC gap 
when connecting two Hamiltonians of different topology 
(winding number) in a continuous manner. To verify this, 
we compute the Floquet-Bloch bands of the Hamiltonian 

( 22 ) 

where is the Hamiltonian describing the sur¬ 

face of the topological insulator in the presence of a 
circularly polarized laser of frequency Q whose helicity 
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FIG. 6 . (color online), (a) Solution map of Eq. 0 [z.e., ln(l + |detQ|)] for the case of two laser beams with different 
frequencies. The used parameters are: = $1,^2 = 2Q, 'y/hO, = 0.2 and (p = 0. (b) Zone center gap closing as a function 

of the relative strength a between the lasers. In the limit case a = —1 (+1) the system is being irradiated with a single 
laser of frequency (Q 2 ) and helicity r = +1 (—1). The parameters coincide with those in panel (a). Inset: Scheme of the 
considered Floquet Hamiltonian indicating the different replicas and their couplings given by the lasers. The lines involve a 
photon absorption/emission of energy hQ (red lines) or 2hQ (blue lines) associated to each one of the lasers. 


is determined by r = ±1. Here, the constant a sim¬ 
ply parametrizes the continuous change of H{t) from 
i^+(f2i,t) to i7_(f22,t). In the case described in the 
main manuscript, this parameter would depend on the 
coordinate perpendicular to the interface and represents 
the relative strength between the two laser beams. Here, 
however, it is an interpolation parameter and thus H{t) 
can be regarded as a local Hamiltonian describing the 
interface region. By taking a particular symmetry point 
(k = 0), we find that when a goes from -1 to 1, hence 
connecting the two topological phases, the ZC gap neces¬ 
sarily closes as indicated in Fig.j^b). This occurs for the 
value of a indicated by the arrow when the replicas from 
m = —2 to m = 2 are retained (the specific value of a 
where the gap closes shows little change when consider¬ 
ing more replicas). Although not shown, the ZC gap also 
closes in other k-points as a is swept. Thus, the gap clos¬ 
ing at certain points in the /c-space as connecting the two 
topological distinct phases ensures bulk-boundary corre¬ 
spondence to hold in this system. Very importantly, TRS 
is always broken regardless the value of a. 
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